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of Argentina: temporal and regional trends 

Ernesto Shimizu^ Analia Macias^, Fernando PaolicchP, Gabriel Magnano^, 
Laura Zapata^, Analia Fernandez^, Ana Canal"*, Sergio Garbaccio^, 
Angel CataldiV^ Karina Caimi"", Martm Zumarraga'' 

'Bacteriology Laboratory, School of Agricultural Sciences, National University of Mar del Plata, Buenos Aires, Argentina 
^School of Agronomy and Veterinary, National University of Rfo Cuarto, Cordoba, Argentina ^General Direction of Animal Production, 
Secretary of Production, Entre Rfos, Argentina ^Department of Basic Pathology, School of Veterinary Sciences, National University of Literal, 
Santa Fe, Argentina 'Pathobiology Institute ''Biotechnology Institute, Centro de Investigacion en Ciencias Veterinarias y Agronomicas, 
Instituto Nacional de Tecnologfa Agropecuaria, Castelar, Buenos Aires, Argentina 

Mycobacterium bovis is the causative agent of bovine tuberculosis (TB), a disease that affects approximately 
5% of Argentinean cattle. Among the molecular methods for genotyping, the most convenient are spoligotyping and 
variable number of tandem repeats (VNTR). A total of 378 samples from bovines with visible lesions consistent with 
TB were collected at slaughterhouses in three provinces, yielding 265 M. bovis spoligotyped isolates, which were 
distributed into 35 spoligotypes. In addition, 197 isolates were also typed by the VNTR method and 54 combined 
VNTR types were detected. There were 24 clusters and 27 orphan types. When both typing methods were combined, 
98 spoligotypes and VNTR types were observed with 27 clusters and 71 orphan types. By performing a meta-analysis 
with previous spoligotyping results, we identified regional and temporal trends in the population structure ofM. 
bovis. For SB0140, the most predominant spoligotype in Argentina, the prevalence percentage remained high during 
different periods, varying from 25.5-57.8% (1994-2011). By contrast, the second and third most prevalent spoligo- 
types exhibited important fluctuations. This study shows that there has been an expansion in ancestral lineages as 
demonstrated by spoligotyping. However, exact tandem repeat typing suggests dynamic changes in the clonal popu- 
lation of this microorganism. 
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Bovine tuberculosis (BTB) is an infectious disease 
primarily caused by Mycobacterium bovis, which affects 
a wide range of mammals, including humans (Amanfu 

2006) . M. bovis is a member of the Mycobacterium tu- 
berculosis complex (MTBC), which also includes M. tu- 
berculosis, Mycobacterium africanum, Mycobacterium 
canetti, Mycobacterium microti, Mycobacterium caprae 
and Mycobacterium pinnipedii (Cataldi & Romano 

2007) and the recently described species Mycobacte- 
rium mungi (Alexander et al. 2010) and Mycobacterium 
orygis (van Ingen et al. 2012). 

In Argentina, the rate of condemnation of carcasses 
exhibiting lesions consistent with BTB dropped from 
6.7% to 0.7% between 1969-2011 (de Kantor et al. 2012). 
The last national screening for purified protein derivative 
(PPD) reactors was performed in 1972 (Nader & Husberg 
1988). Therefore, the current prevalence of PPD reactors 
is unknown, but it may be higher than 0.7%i and even 
higher in dairy cattle. The control strategies are based on 
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the caudal fold skin test, which is used primarily in dairy 
herds and slaughter surveillance. M. bovis has been iso- 
lated from tuberculin-positive bovines after a post mor- 
tem inspection of reactors in the slaughterhouse. 

Although there are no reports of a phenotypic method 
to distinguish between different M bovis isolates, geno- 
typing methods have achieved some distinction. Among 
the molecular methods available for genotyping, the 
most convenient are spoligotyping (a reverse line probe 
assay) and genotyping by variable number of tandem re- 
peats (VNTR). 

At present, the methods for genotyping M. bovis and 
other MTBC species include restriction fragment length 
polymorphism (Otal et al. 1991, Collins et al. 1993) 
and VNTR (Haddad et al. 2001). However, until now, 
spoligotyping has been the best option for large-scale 
screening studies on the distribution of M bovis strains 
(Kamerbeek et al. 1997). The theories of molecular fin- 
gerprinting suggest that epidemiologically related iso- 
lates will have similar fingerprints, which are expected 
to be different from epidemiologically unrelated strains 
(Maslow et al. 1993). Desirable characteristics for a typ- 
ing method are related to the stability of the strain and 
the diversity within the species (Skuce & Neill 2001). 
Although spoligotyping may fulfil these characteristics 
in some cases, other typing methods such as VNTR must 
be performed to discriminate among isolates, especially 
in settings where one spoligotype is largely predominant 
or has been used to investigate a recent transmission. 
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Previous studies have demonstrated a wide diversity 
of M bovis spoligotypes for cattle populations in Latin 
American countries. For example, different M. bovis 
spoligotypes have been found in different regions of Ar- 
gentina, Brazil, Uruguay and Mexico (Zumarraga et al. 
1999, Milian-Suazo et al. 2002, Rodriguez et al. 2004). 
The distribution of circulating M. bovis spoligotypes in 
five Latin American countries has been analysed recent- 
ly (Zumarraga et al. 2012). 

The aim of this work was to assess the genetic diver- 
sity of M. bovis in Argentina and to envisage temporal 
and spatial trends in the clonal population of this patho- 
gen in the context of active transmissions in cattle. 

MATERIALS AND METHODS 

Animals - The bovines used in this study were exam- 
ined in slaughterhouses in three provinces of Argentina, 
namely Buenos Aires, Cordoba and Santa Fe. The ani- 
mals were classified according to age and weight as fol- 
lows: calves 10% (< 12 months, < 220 kg), young steers 
5% (12-18 m < 350 kg), steers 14% (> 18 m, > 350 kg), 
heifers 8.5% (12-30 m, 260 kg), cows 60.5% (> 30 m, > 
350 kg) and bulls 0.5% (> 24 m). 

Bacterial cultures - Tissue samples collected from 
378 bovines with lesions classified as tuberculous were 
detected during regular post mortem inspections in na- 
tional abattoirs. 

A convenience, non-probabilistic sampling was per- 
formed at each slaughterhouse with the assistance of the 
local animal health service staff (National Health Serv- 
ice and Food Quality). For bacteriological analysis, the 
samples were placed in sterile receptacles and stored at 
-20°C until processing. Petroff's method was used to 
decontaminate the samples for mycobacterial isolates. 
Each sample was inoculated on Stonebrink and Lowen- 
stein- Jensen media at 37°C for 60 days. Each culture was 
performed in duplicate. During this period, the cultures 
were examined every seven days and when growth was 
visible, smears were prepared and stained with Ziehl- 
Neelsen to observe acid-fast bacilli. The bacteriological 
identification of the isolates was performed on the basis 
of the culture media, incubation temperature, growing 
time, colony morphology and Ziehl-Neelsen staining. 

Isolates - A total of 261 M. bovis isolates were includ- 
ed in this study. The isolates originated from bovines at 
slaughterhouses in Buenos Aires (n = 67), Cordoba (n = 
112) and Santa Fe (n = 82). These isolates were obtained 
between 2008-2011. These three Argentinean provinces 
are inhabited by 60% of the cattle in the country. 

Spoligotyping assays and analysis - Spoligotyping 
was performed according to Kamerbeek et al. (1997). 
Bovine tissue samples were decontaminated by using 
Petroff's method and cultured in Lowenstein-Jensen 
and Stonebrink media at 37''C for 60 days. A loopful 
of colonies was transferred into a microcentrifuge tube 
containing 250 (xL of distilled water and heated at 96°C 
for 45 min. The samples were then centrifuged at 12,000 
rpm for 10 min and 5 |iL of the supernatant was used 
for polymerase chain reaction (PCR) to amplify the DR 



region. Spoligotyping was performed with a spoligotyp- 
ing kit (Ocimun Biosolutions, Hyderabad, India). M. tu- 
berculosis H37Rv (ATCC 27294) and M. bovis Bacillus 
Calmette-Guerin (BCG) (ATCC 27289) were included 
as reference strains in each spoligotyping experiment. 
The spoligotypes were collected in binary format in an 
Excel database and scanned images of the films were 
analysed by using BioNumerics (v.3.5. Applied Maths, 
Sint-Martens-Latem, Belgium). Spoligotypes were com- 
pared with the M. bovis spoligotypes contained in the 
Biotechnology Institute of National Institute of Agricul- 
tural Tecnology database, Argentina, and in the mbovis. 
org database from the Veterinary Laboratories Agency 
(VLA) (servers hosted at the University of Sussex in the 
United Kingdom). The spoligotypes were named accord- 
ing to VLA code (SB). We defined the following criteria 
for the purpose of this study: a unique spoligotype or 
VNTR type was observed only once within a province; 
an orphan spoligotype was observed only once during 
the study and an exclusive spoligotype was clustered or 
not clustered, in only one province. 

VNTR typing - The analysis of the six VNTR loci 
originally identified by Frothingham and Meeker- 
O'Connell (Frothingham & Meeker-O'Connell 1998) 
was limited to the exact tandem repeat (ETR) of ETR-A 
to -D loci, as previously performed by Smith et al. (2003). 
The VNTR genotype of a strain, which represents the 
number of repeat elements at each locus, is presented as 
a series of four integers between 1-12 that are separated 
by hyphens. When the variant of an integer was identi- 
fied, it was marked by an asterisk (*). 

Multiplex PCRs were used to combine primer pairs 
ETR-A/B and ETR-C/D. The PCR mix was prepared in 
96 -well plates using the Hot Start Master Mix Kit (Qia- 
gen, Germany). Five nanograms of DNA were added to a 
final volume of 20 (xL containing 0.4 (xM of each primer. 
For each multiplex mixture, one primer of each oligonu- 
cleotide pair was tagged with a different fluorescent dye 
(Table I). The thermocycler programs for both multiplex 
reactions were identical. The PCRs were performed with 
an initial denaturation of 15 min at 94°C followed by 35 
cycles of denaturation at 94°C for 30 s, annealing at 68°C 
for 1 min and extension at 72''C for 2 min, followed by a 
final extension at 72°C for 10 min. 

Statistical analysis - The discriminatory index (D) 
described by Hunter and Gaston and expressed by the 
Simpson formula was calculated to determine the dis- 
criminatory power of spoligotyping in each province 
(insilico.ehu.es) (Rodriguez et al. 2010). This index is 
given by the following equation (N: total number of 
strains in the sample population; s: total number of types 
described; xj: number of strains belonging to the j* type) 
(Hunter & Gaston 1988): 



This index is based on the probability that two unrelat- 
ed strains sampled from the test population will be placed 
into different typing groups (Hunter & Gaston 1988). 



238 Mem Inst Oswaldo Cruz, Rio de Janeiro, Vol. 709(2), April 201 4 



TABLE I 

Primers for exact tandem repeat (ETR)-A-D typing 



Primer sequence (5-3) 

Variable number of tandem repeats locus (labelling) 

ETR-A AAATCGGTCCCATCACCTTCTTA-(FAM) 

CGAAGCCTGGGGTGCCCGCGATTT 
ETR-B GCGAACACCAGGACAGCATCAT 

(JOE)-GGCATGCCGGTGATCGAGTGG 
ETR-C GTGAGTCGCTGCAGAACCTGCAG-(HEX) 

GGCGTCTTGACCTCCACGAGTG 
ETR-D CAGGTCACAACGAGAGGAAGAGC-(FAM) 

GCGGATCGGCCAGCGACTCCTC 



RESULTS 

Geographical distribution of isolates and spoligo- 
types - Valid spoligotypes were obtained for 261 isolates 
across 32 patterns. Of these spoligotypes, 13 patterns 
were orphans and 248 (94.6%) isolates were grouped in 
18 clusters of more than two isolates (Table II), with a 
discrimination index of 0.80. The spoligotype SB0140 
represented 39.8% of the isolates. The isolates from Cor- 
doba specifically came from slaughterhouses located in 
the Juarez Celman department. The numbers of spoli- 
gotypes detected in Buenos Aires, Cordoba and Santa 
Fe were nine, 21 and 19, respectively. A map showing 
the distribution of the four major spoligotypes is shown 
in Fig. lA. These data can be compared to the cattle 
density per department (Fig. IB) and it should be noted 
that in general samples were taken from departments 
with higher numbers of cattle. In Buenos Aires, SB0140 
comprised 28.4% (n = 19) of the isolates, SB0145, 35.8% 
(n = 24), SB0120, 16.4% (n = 11) and SB1055, 10.4% (n 
= 7). In Cordoba, SB0140 comprised 45.5% of the iso- 
lates (n = 51), SB0145, 8.9% (n = 10) and SB0153 and 
SB0484, 8% (n = 9 each). In Santa Fe, which has a strong 
dairy cattle industry, SB0140 comprised 41.4% (n = 34), 
SB0130, 14.6% (n = 12) and SB0120, 10.9% (n = 9). Bue- 
nos Aires had the lowest number of unique spoligotypes 
(n = 4) followed by Cordoba (n = 11) and Santa Fe (10). 
Regarding orphan spoligotypes, Buenos Aires presented 
only one orphan spoligotype (xxl38), Cordoba had sev- 
en (SB0996, 1043, 1068, 1776, xxl26, xxl31 and xxl34) 
and Santa Fe had six (SB0678, 1066, 2117, 2165, 2166 
and xxl39). The number of major spoligotypes, the re- 
maining clustered spoligotypes and the orphan types are 
shown in Fig. 2. 

The exclusive spoligotypes were composed of the 
orphans in addition to four clustered spoligotypes, with 
one from Buenos Aires (SB 1055, n = 7), three from 
Cordoba (SB0269, n = 2, SB0273, n = 4 and SB0520, 
n = 4) and one from Santa Fe (SB1790, n = 2). The dis- 
tribution of spoligotypes in each province is shown in 
Table II. Four spoligotypes (SB0120, SB0131, SB0140 
and SB0145) were shared among the three provinces 
and they made up 65.9%) of the isolates. Six spoligo- 



TABLE II 

Spoligotypes distribution, number of patterns per province 



Spoligotype 


Total 


Buenos 




Santa 


(SB code) 


of isolates 


Aires 


Cordoba 


Fe 


SB0120 


24 


11 


4 


9 


SB0130 


18 




6 


12 


SB0131 


8 


2 


2 


4 


SB0140 


104 


19 


51 


34 


SB0145 


36 


24 


10 


2 


SB0153 


11 




9 


2 


SB0267 


2 


- 


1 


1 


SB0269 


2 


- 


2 


- 


SB0273 


4 




4 




SB0484 


13 




9 


4 


SB0520 


4 




4 




SB0678 


1 






1 


SB0820 


2 


1 


1 




SB0980 


4 




1 


3 


SB0996 


1 




1 




SB1033 


2 


1 




1 


SB1043 


1 




1 




SB1055 


7 


7 






SB1066 


1 






1 


SB1068 


1 




1 




SB1776 


1 




1 




SB1787 


2 


1 




1 


SB1790 


2 






2 


SB1799 


2 




1 


1 


SB2117 


1 






1 


SB2165 


1 






1 


SB2166 


1 






1 


SBxxl26 


1 




1 




SBxxBl 


1 




1 




SBxxl34 


1 




1 




SBxxl38 


1 


1 






SBxxl39 


1 






1 


Total 


261 


67 


112 


82 
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types were present in Cordoba and Santa Fe (SB0130, 
SB0153, SB0267, SB0484, SB0980 and SB1799), which 
represented 26.8% of the isolates from these provinces. 
Another spoligotype was detected in Cordoba and Bue- 
nos Aires (SB0820) (with 1.1% of the isolates from these 
2 provinces) and two spoligotypes were in Buenos Aires 
and Santa Fe (SB1033 and SB1787) (with 2.7% of the iso- 
lates from these 2 provinces). 

With regards to the distribution of spoligotypes per 
animal category, we observed that in the calf category 
(n = 27), SB0145 was predominant, with 19 isolates. The 
animals from this specific category were primarily sam- 
pled from Buenos Aires (81.5%). There were 13 samples 
in the category for young steers. All of these samples 
came from Cordoba and SB0140 was the predominant 
spoligotype (n = 5) in this category. Samples from 35 
animals were included in the steer category (65.7%) 
from Cordoba and 31.4%o from Santa Fe). Out of these 
samples, 16 belonged to SB0140 and five to SB0145. In 
the heifer category, most animals came from Cordoba 
(83.3%). Of the 24 animals from this particular category, 
nine belonged to SB0140 and six to SB0145. In the cow 
category (n = 158) (Buenos Aires, n = 40; Cordoba, n 
= 49 and Santa Fe, n = 69), the following spoligotypes 
were described: SB0120 (19), SB0130 (13), SB0140 (68), 
SB0153 (5), SB0484 (11) and SB1055 (7). Only one ani- 
mal, which came from Cordoba, was included in the bull 
category and it exhibited spoligotype SB0140. In Santa 
Fe, there were 11 isolates from the steer category and 
spoligotype SB0140 included 45.5% of the isolates. From 
cows, we obtained 69 isolates and the spoligotypes were 
distributed as follow: SB0140 (28), SB0130 (12), SB0120 
(7) and SB0484 (4). In Cordoba, SB0140 predominated 
in all categories [young steers (38.5%)), steers (47.8%o), 
heifers (45%) and cows (49%)]. 

VNTR typing - A total of 197 isolates were also typed 
by VNTR (Buenos Aires, n = 64, Cordoba, n = 69 and 
Santa Fe, n = 64). There were 55 combined VNTR types 
that were determined by using primers directed to ETR 



A, B, C and D (Table III). The discrimination index was 
0.92. There were 25 clusters (comprising 168 isolates) 
and 29 orphan types. The predominant cluster was 7-5- 
5-4* with 48 isolates present in the three provinces. The 
second most prevalent cluster was 6-5-5-4* with 20 iso- 
lates present in the three provinces, followed by 6-5-5-1 
with 12 isolates from Santa Fe. 

Subtyping of spoligotypes by VNTR - The highly 
predominant spoligotype SB0140 was divided in 26 
VNTR types, where 7-5-5-4* included 15 of the 64 
isolates (data not shown). The 34 isolates belonging to 
SB0145 were divided into eight VNTR types. In turn, 
the third most frequent spoligotype, namely SB0120 
(with 21 isolates), was subdivided into 11 VNTR types. 
With regards to the major VNTR types, the predomi- 
nant 7-5-5-4* with 48 isolates was split into 14 differ- 



Total (numb«i) Buenos Aires (number) 




Santa Fe (number) 

Cordoba (number) 




Fig. 2: pie chart representing tlie major spoligotypes, tlie rest of clus- 
tered spoligotypes and orphan types in the three provinces and in 
each province. 
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Fig. lA: distribution per department of the four major spoligotypes in the three provinces (1 dot = 1 case). Department are the secondary admin- 
istrative and political division of provinces in Argentina; B: number of cattle per department (source: rian.inta.gov.ar/ganaderia/stock2011). 
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TABLE III 

Variable number of tandem repeats (VNTR) 
types distribution, number of patterns per province 



VNTR 


Total 


Buenos 




Santa 


types 


of isolates 


Aires 


Cordoba 


Fe 


2-3-5-3 


1 


1 


- 


- 


2,5*-5-4,* 


1 


1 


- 


- 


3-10-5-4* 


1 


- 


- 


1 


3-3-5-1 


5 


5 


- 


- 


3-5-1-4* 


1 


1 


- 


- 


3-5-5-1 


1 


1 


- 


- 


3-5-5-4* 


1 


3 


1 


- 


4-4-5-1 


1 


1 


- 


- 


4-5-5-1 


3 


2 


- 


1 


4-5-5-4* 


3 


- 


2 


1 


5-4-5-1 


1 


- 


- 


1 


5-5-1-4* 


2 


- 


- 


2 


5-5-5-4* 


3 


- 


2 


1 


5-8-4-4* 


1 


1 


- 


- 


6-1-5-4* 


1 


- 


1 


- 


6-2-1-4* 


1 


- 


1 


- 


6-3-5-1 


3 


- 


- 


3 


6-3-5-4* 


1 


- 


- 


1 


6-3*-5-4* 


2 


- 


1 


1 


6-4-5-4* 


1 


- 


1 


- 


6-5-4-4* 


1 


- 


1 


- 


6-5-5-1 


12 


- 


- 


12 


6-5-5-3 


2 


1 


- 


1 


6-5-5-4* 


20 


4 


12 


4 


6-9-4-4* 


1 


- 


- 


1 


6-9-5-1 


2 


- 


- 


2 


6-9-5-3 


1 


- 


- 


1 


6-9-5-4* 


6 


4 


- 


2 


6-9*-5-4* 


1 


1 


- 


- 


7-10-4-1 


1 


- 


1 


- 


7-10-4-3 


1 


- 


1 


- 


7-10-4-4* 


1 


- 


- 


1 


7-10-5-1 


6 


1 


- 


5 


7-10-5-3 


1 


- 


- 


1 


7-10-5-4* 


2 


1 


- 


1 


7-3-5-4* 


1 


- 


- 


1 


7-3*-5-4* 


1 


1 


- 


- 


7-4-4-1 


1 


- 


- 


1 


7-4-5-1 


5 


2 


- 


3 


7-4-5-4* 


6 


4 


1 


1 


7-5-1-4* 


2 


- 


1 


1 


7-5-2*-4 


2 


2 


- 


- 


7-5-3-1 


1 




- 


1 


7-5-4-4* 


1 


1 






7-5-5-1 


5 




4 


1 


7-5-5-3 


9 




2 


1 


7-5-5-4* 


48 


10 


29 


9 


7-5-6-4* 


2 




2 




7-5-7*-4* 


1 




1 




7-6-4-1 


1 






1 


7-6-5-4* 


4 




4 




7-7-5-4* 


1 




1 




7-9-5-1 


1 






1 


7-9*-5-4* 


10 


10 






Total 


197 


64 


69 


64 



ent spoligotypes, which were primarily distributed into 
the following three spoligotypes: SB0153, SB0145 and 
SB0140. The second most prevalent VNTR was 6-5- 
5-4*, with 20 isolates composing spoligotype SB0140 
in 14/20 isolates and this VNTR was split into six dif- 
ferent spoligotypes. Finally, 6-5-5-1 comprised SB0140 
and SB0130 with 12 isolates each. 

Combined discrimination byspoligotypingand VNTR 
typing - When both typing methods were combined, we 
detected 96 spoligotypes and VNTR types (Table IV) 
with 30 clusters (comprising 131 isolates) and 66 orphan 
types and a discrimination index of 0.98. SB0140-7-5- 
5-4* was the predominant cluster with 15 isolates dis- 
tributed in South Cordoba, Santa Fe and Buenos Aires, 
followed by SB0140-6-5-5-4* with 14 isolates distrib- 
uted in the three provinces. Then, SB0145-7-5-5-4* was 
present in Cordoba and Buenos Aires and SB0145- 7-9*- 
5-4* was observed only in south Buenos Aires with 10 
isolates each. 

In addition to the orphan types, there were 18 spo- 
ligo-VNTR types that were exclusives (6 of each prov- 
ince), which were clustered in 60 isolates (Buenos Aires, 
25, Cordoba, 15 and Santa Fe, 20). 

Temporal trends - A study of temporal trends was 
performed with a meta-analysis of spoligotypes. We 
used data from 1994-2011 that were recorded in our 
genotyping database. The analysis was performed by di- 
viding the data into periods of four years (except for the 
1st period that includes only the years 1994 and 1995). 
The predominant spoligotype was SB0140 and its preva- 
lence percentage remained high throughout the different 
periods, varying from 25.5-57.8% (Fig. 3). SB0130 was 
the second most prevalent in the database registers and 
experienced an important fluctuation and it reached a 
peak of 65.2% in 2000-2003. Similarly, SB0145 also dis- 
played fluctuations, although they were at lower preva- 
lence levels. By contrast, other minor spoligotypes had 
low representation, but they were constant or sporadic. 
For instance, SB0131 had low representation and SB0484 
was sporadically represented. 

When temporal trends are taken into account by 
province, the most predominant spoligotype in Buenos 
Aires was SB0140, which is in line with the rest of the 
country. However, the prevalence percentage highly dif- 
fers throughout the different time periods (Fig. 3), from 
100-16%). SB0145, the second most predominant type in 
Buenos Aires, increased sharply during the last study 
period. Additionally, important fluctuations are ob- 
served in the periodical frequency of SB0130. In some 
cases, the sampling number was too low to validate the 
analysis (data not shown). This challenge was also the 
case for the analysis in Cordoba, where the number of 
collected isolates was too low during some periods (data 
not shown). In Santa Fe, SB0140 was the predominant 
spoligotype with a percentage above 37%o during all 
periods, except for 2000-2003. The second most pre- 
dominant spoligotype was SB0130, which varied from 
0-93.75%. Finally, SB0120 rose from 0-10.2% in 2008- 
2011 (data not shown). 
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TABLE IV 

Distribution of combined spoligotyping-variable number of tandem repeats (VNTR) types, number of patterns per province 







Total 


Buenos 




Sant 


Spoligotype 


VNTR 


of isolates 


Aires 


Cordoba 


Fe 


SB0120 


2- 3- 5- 3 


1 


1 


- 


- 


SB0120 


3- 5- 5- 1 


1 


1 


- 


- 


SB0120 


5- 5- 1-4* 


2 


- 


- 


2 


SB0120 


5-5-5-4* 


1 


- 


1 


- 


SB0120 


7- 10- 5- 4* 


1 


1 


- 


- 


SB0120 


7- 3- 5- 4* 


1 


- 


- 


1 


SB0120 


7- 4- 5- 1 


3 


2 


- 


1 


SB0120 


7- 4- 5- 4* 


5 


4 


- 


1 


SB0120 


7-5-5-3 


1 


- 


1 


- 


SB0120 


7- 5- 5- 4* 


4 


2 


1 


1 


SB0120 


7- 6- 5- 4* 


1 


- 


1 


- 


SB0130 


6- 2- 1- 4* 


1 


- 


1 


- 


SB0130 


6-3-5- 1 


2 


- 


- 


2 


SB0130 


6- 3*- 5- 4* 


2 


- 


1 


1 


SB0130 


6-5-5- 1 


6 


- 


- 


6 


SB0130 


6-5-5-4* 


2 


- 


2 


- 


SB0130 


6-9-5-4* 


1 


- 


- 


1 


SB0130 


7- 5- 5- 4* 


3 


- 


2 


1 


SB0131 


3- 10- 5-4* 




- 


- 


1 


SB0131 


4-4-5- 1 




1 


- 


- 


SB0131 


7- 4- 5- 4* 




- 


1 


- 


SB0131 


7- 5- 3- 1 




- 


- 


1 


SB0131 


7- 5- 5- 1 




- 


- 


1 


SB0131 


7- 5- 5- 4* 




1 


- 


- 


SB0140 


2*- 5- 5- 4* 




1 


- 


- 


SB0140 


3-5-5-4* 


2 


1 




- 


SB0140 


4-5-5- 1 


2 


2 




- 


SB0140 


4-5-5-4* 


2 


- 




1 


SB0140 


5-4-5- 1 




- 




1 


SB0140 


5-5-5-4* 




- 




- 


SB0140 


5- 8- 4- 4* 




1 




- 


SB0140 


6- 1-5-4* 




- 




- 


SB0140 


6-3-5- 1 




- 




1 


SB0140 


6-5-4- 4* 




- 




- 


SB0140 


6-5-5- 1 




- 


- 


6 


SB0140 


6-5-5-3 




- 


- 


1 


SB0140 


6-5-5-4* 


14 


4 


6 


4 


SB0140 


6- 9-4- 4* 




- 


- 


1 


SB0140 


6-9-5-3 




- 


- 


1 


SB0140 


6-9-5-4* 




4 


- 


- 


SB0140 


6- 9*- 5- 4* 




1 






SB0140 


7- 10- 4- 4* 








1 


SB0140 


7- 10- 5- 4* 








1 


SB0140 


7- 3*- 5- 4* 




1 






SB0140 


7- 4- 4- 1 








1 


SB0140 


7- 5- 4- 4* 




1 






SB0140 


7- 5- 5- 3 








1 


SB0140 


7- 5- 5- 4* 


15 


1 


10 


4 
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Total 


Buenos 




Santa 


Spoligotype 


VNTR 


of isolates 


Aires 


Cordoba 


Fe 


SB0140 


7- 5- 7*- 4* 


1 


- 


1 


~ 


SB0140 


7- 6- 4- 1 


1 


- 




1 


SB0145 


3- 5- 1- 4* 


1 


1 


~ 


~ 


SB0145 


6- 5- 5- 4* 


1 


~ 


1 


~ 


SB0145 


7- 4- 5- 1 


1 


~ 




1 


SB0145 


7- 5- 2*- 4 


2 


2 


~ 




SB0145 


7- 5- 5- 3 


7 


6 


1 




SB0145 


7- 5- 5- 4* 


10 


5 


5 




SB0145 


7- 6- 5- 4* 


2 


~ 


2 




SB0145 


7- 9*- 5- 4* 


10 


10 


- 


~ 


SB0153 


6- 3- 5- 4* 


1 


~ 


- 


1 


SB0153 


6- 5- 5- 4* 


1 




1 




SB0153 


7- 5- 1- 4* 


1 




1 




SB0153 


7- 5- 5- 4* 


4 




3 


1 


SB0153 


7- 5- 6- 4* 


2 




2 


~ 


SB0153 


7- 6- 5- 4* 


1 




1 




SB0153 


7- 7- 5- 4* 


1 




1 


~ 


SB0267 


7- 5- 5- 4* 


1 




- 


1 


SB0269 


6- 4- 5- 4* 


1 




1 




SB0273 


7- 10- 4- 1 


1 




1 




SB0273 


7- 10- 4- 3 


1 




1 




SB0273 


7- 5- 5- 4* 


2 




2 


~ 


SB0484 


4- 5- 5- 1 


1 


- 


- 


1 


SB0484 


7- 10- 5- 1 


2 


- 


- 


2 


SB0484 


7- 5- 5- 1 


4 


- 


4 


- 


SB0484 


7- 5- 5- 4* 


1 


- 


1 


- 


SB0520 


7-5-5-4* 


3 




3 




SB0820 


7- 10- 5- 1 


1 


1 


~ 


~ 


SB0980 


5- 5- 5-4* 


1 






1 


SB0980 


6- 9- 5- 1 


1 






1 


SB0980 


7- 9- 5- 1 


1 






1 


SB0996 


6- 5- 5- 4* 


1 




1 


~ 


SB1033 


6- 5- 5- 3 


1 


1 


~ 


~ 


SB1033 


6- 9- 5- 1 


1 






1 


SB1043 


7- 5- 5- 4* 


1 




1 


~ 


SB1055 


3- 3- 5- 1 


5 


5 






SB1055 


3- 5- 5-4* 


2 


2 


~ 




SB1068 


4- 5- 5- 4* 


1 


~ 


1 




SB1787 


7-4-5- 1 


1 


1 




~ 


SB1787 


7- 5- 1- 4* 


1 






1 


SB1790 


7- 10- 5- 1 


2 






2 


SB1799 


7- 5- 5- 4* 


1 




1 




SB2117 


7- 10- 5- 1 


1 






1 


SB2165 


6- 9- 5- 4* 


1 




. 


1 


SB2166 


7- 10- 5- 3 


1 






1 


Spol26 


6-5-5-4* 


1 




1 




Spol38 


7-5-5-4* 


1 


1 






Spol39 


7-5-5-4* 


1 






1 


Total 




197 


64 


69 


64 
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■ 2008-201 1 

■ 2007-2004 

□ 2003-2000 

□ 1999-1996 

■ 1995-1994 



Fig. 3: percentage of incidence of major spoligotypes in the three 
studied provinces along different sampling periods. 



DISCUSSION 

Among the 32 spoligotypes detected, 10 have not been 
previously included or described in Argentina and are 
not included in the databases from mbovis.org; therefore, 
they were considered to be new spoligotypes. Of these 10 
spoligotypes, three formed clusters of two isolates each, 
that is, SB1790 of Santa Fe, SB1787 of Santa Fe and Bue- 
nos Aires and SB1799 of Santa Fe and Cordoba. The oth- 
er new spoligotypes were also orphans, namely SBxxl38 
from Buenos Aires, SBxxlBl and 1776 from Cordoba 
and SB2165, 2166, xxl26 and xxl39 from Santa Fe. Santa 
Fe had the majority of new spoligotypes, most likely be- 
cause of the active TB eradication program in that prov- 
ince. This program, which was implemented across the 
whole country, consists of the sacrifice of infected bo- 
vines and the inspection of carcasses in slaughterhouses. 
In this way, when a new M bovis strain is introduced 
into a herd, the infected animals are detected early and 
sacrificed, so a new M. bovis genotype can be detected 
before it can spread among the other animals. Other au- 
thors have proposed that if a human or bovine TB control 
program is active and long-lasting, then a reduction in 
predominant clones should be observed as transmission 
is halted (van Soohngen et al. 1991, Durr et al. 2000). 

Early in the application of VNTR elements to M bo- 
vis, the method's discriminatory power was observed to 
be higher than that of spoligotyping (Roring et al. 2004). 
Specifically, ETR A and ETR B exhibited an important 
allelic variation in M. bovis (Romero et al. 2011). The 
ETR A-C elements were previously shown to have a dis- 
crimination index higher than 0.4 (Martinez et al. 2008, 
Boniotti et al. 2009, Price-Carter et al. 2011). However, 
the individual discrimination indexes in M. bovis for 
ETR A-D may be highly variable according to the geo- 
graphic setting (Sahraoui et al. 2009, Sun et al. 2012). 
This discrimination index increases when ETRs A-D 
are evaluated together. Additionally, the combination of 
VNTR typing and spoligotyping increases the discrimi- 
nation index above that of VNTR alone (Martinez et al. 
2008). In contrast with the ETR- VNTR elements, mini- 
mum inhibitory concentration (MIRU)-VNTR (with 
the exception of MIRU26) yields a low discriminatory 
power for M. bovis (Rodriguez et al. 2004, Parreiras et 
al. 2012, Rocha et al. 2013). 



One important aspect to evaluate is the discrimination 
conferred by VNTR typing in predominant spoligotypes. 
We observed that the predominant SB0140 spoligotype 
was split into 26 different ETR A-D types. This result is 
in accordance with the findings of McLernon et al. (2010) 
in Ireland, where SB0140 is also the predominant spoli- 
gotype and was highly subtyped by VNTR typing. Bon- 
iotti et al. (2009) observed that 408 isolates of SB0120, 
which is the most predominant spoligotype in Italy, were 
split into 64 different patterns (Boniotti et al. 2009). 

The presence of combined cluster patterns is an index 
for the active transmission of TB. There were 66.5% of 
the isolates distributed into 30 different spoligotyping- 
VNTR combined patterns. Smith et al. (2003) suggested 
that the clonal expansion of a genotype is caused either by 
the spread of a favourable mutation or by the invasion of a 
geographic region by a few genotypes (Smith et al. 2003). 

In this work, we detected 18 combined and clustered 
spoligotyping-VNTR patterns (6 in each province) that 
were exclusives in the sampled provinces. These com- 
bined patterns represented 30.5% of the analysed isolates. 
Together with the 66 orphans, these patterns yielded 84 
combined and exclusive patterns. This finding suggests 
a pattern association with a particular geographical re- 
gion. However, this putative association must be demon- 
strated with more isolates in the future; the association 
was supported by a chi-square test because of the low 
number of isolates (data not shown). 

Our database consists of M. bovis molecular typ- 
ing data from 1994 until the present. When we divided 
this 17-year span into regular periods, we observed that 
although the incidence of different spoligotypes fluctu- 
ated, the most predominant spoligotypes were unaltered. 
To the best of our knowledge, only one previous publica- 
tion on the evolution of the M. bovis population structure 
has been reported (Haddad et al. 2001). 

The predominance of major spoligotypes in cattle 
from the primary cattle farms in these Argentine prov- 
inces suggests that there is still an active transmission 
of the pathogen, in spite of the reduced incidence of 
animals with visible lesions in slaughterhouses. Spoli- 
gotypes are more stable than VNTR types. In fact, the 
BCGs obtained in 1921 from which the substrains origi- 
nated cannot be differentiated by spoligotyping (de la 
Salmoniere et al. 1997, Bauer et al. 1999) and VNTR 
subdivides BCG (Magdalena et al. 1998, Supply et al. 
2000). As a consequence, the predominant spoligotypes 
may ancient, as suggested by the subdivision of the main 
spoligotypes by VNTR typing. 
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